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The ground state of the uniform antiferromagnetic spin-1/2 Heisenberg chain can be viewed as 
strongly fluctuating liquid of valence bonds, while in disordered chains these bonds lock into random 
singlet states on long length scales. We show that this phenomenon can be studied numerically, even 
in the case of weak disorder, by calculating the mean value of the number of valence bonds leaving 
a block of L contiguous spins (the valence bond entanglement entropy) as well as the fluctuations 
in this number. Results are obtained for a class of models which include the Heisenberg chains as 
special cases. 



The set of valence-bond states — states in which lo- 
cahzed spin-1/2 particles are correlated in singlet pairs 
said to be connected by valence bonds — provides a use- 
ful basis for visualizing singlet ground states of quantum 
spin systems. For example, the ground state of the uni- 
form one-dimensional nearest-neighbor spin-1/2 antifer- 
romagnetic (AFM) Heisenberg model (the prototypical 
spin-liquid state [T]) can be viewed as a strongly fluc- 
tuating liquid of valence bonds with a power-law length 
distribution. This intuitive picture reflects the long-range 
spin correlations in this state, as well as the existence of 
gapless excitations created by breaking long bonds. 

Valence-bond states also play a key role in describing 
the physics of random spin-1/2 AFM Heisenberg chains. 
For these systems, it was shown by Fisher [2j, using a 
real space renormalization group (RSRG) analysis, that 
on long length scales the ground state is described by 
a single valence-bond state known as a random singlet 
state. This single valence-bond state should be viewed as 
a caricature of the true ground state, which will certainly 
exhibit bond fluctuations on short length scales. In fact, 
it is natural to expect that, when measured on these short 
length scales, a fluctuating random singlet state would 
be difficult to distinguish from the uniform Heisenberg 
ground state, particularly in the limit of weak disorder. 

In valence-bond Monte Carlo (VBMC) simulations [5] 
valence-bond states are used to stochastically sample sin- 
glet ground states of quantum spin systems. One of the 
appealing features of VBMC is that if one imagines view- 
ing the sampled valence-bond states over many Monte 
Carlo time steps the resulting "movie" would correspond 
closely to the intuitive resonating valence bond picture 
described above. For random Heisenberg chains (and re- 
lated models) VBMC should therefore provide a useful 
method for directly studying the phenomenon of random 
singlet formation on long length scales, while at the same 
time capturing the short-range fluctuations which will al- 
ways be present. 

With this motivation, we present here results of VBMC 
studies for a class of models which include the uniform 
and random spin-1/2 AFM Heisenberg chains. To de- 
fine these models we first specify their relevant Hilbert 




FIG. 1: (Color online) (a) A non-crossing valence-bond state, 
(b) Action of singlet projection operators n° on a non-crossing 
valence-bond state, (c) A non-crossing valence-bond state 
with a block of L = 5 sites (shaded region) for which n_L, the 
number of bonds leaving the block, is 3. 



spaces. It is well known that the set of non-crossing 
valence-bond states (see Fig. 1(a)) form a complete lin- 
early independent basis spanning the total spin Hilbert 
space of a chain of spin-1/2 particles |3]. We denote the 
singlet projection operator acting on neighboring sites i 
and i -I- 1 by H°, which, for spin-1/2 particles, can be 
expressed as H? = | — 5i • 5*^+1 (taking h = 1). Figure 
1(b) shows two representative examples of H-' acting on a 
non-crossing valence-bond state. For spin-1/2 particles, 
the parameter d appearing in Fig. 1(b) is equal to 2; how- 
ever, in principle, d can take any value, (of course \id^2 
the Hilbert space no longer describes spin-1/2 particles). 
Of particular interest are the cases d = 2 cos where 
fc is a positive integer [S]. For these values of d, when 
k is finite, the non-crossing states are no longer linearly 
independent and the Hilbert space dimensionality of N 
sites can be shown to grow asymptotically as d^ with 
d < 2. The k ^ oo limit then corresponds to the case 
of ordinary spin-1/2 particles with d — 2 for which the 
Hilbert space dimensionality grows as 2^. 

One consequence of the reduced Hilbert space dimen- 
sionality for finite integer k is that it changes the entan- 
glement entropy associated with a valence bond. The 
entanglement entropy of a subsystem ^ of a larger sys- 
tem consisting of parts A and B is defined to be the 
von Neumann entropy, S^^ , of the reduced density ma- 
trix Pa obtained by tracing out the degrees of freedom in 
region B, thus 5*^^ = ~Tr[pAlog2 Pa]- With this defi- 
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nition, an ordinary singlet formed by two spin-1/2 parti- 
cles, with one spin in region A and the other in region B, 
will have S'^^ = 1. However, when d = 2 cos it was 
shown in [B] that if there are M valence bonds connect- 
ing sites in region A with sites in region B, then, in the 
M ^ 1 limit, because the dimensionality of the traced 
out Hilbert space grows as , S^^ ~ A/logg d and the 
entanglement per bond is log2 d. 

The class of Hamiltonians studied here are all charac- 
terized by the parameter d and have the form 

i/ = -^j,no, (1) 

i 

with Ji > 0. For d = 2 these models correspond to 
spin-1/2 AFM Heisenberg chains with Ji equal to the 
exchange energy associated with spins i and i + 1. For 
general d, if the J/s are uniform (J^ = J) the Hamil- 
tonians ([T| can be viewed as 1-1-1 dimensional quantum 
Potts models obtained by taking the asymmetric limit 
of the transfer matrix of the Q-state Potts models with 
Q — d^ [7]. For d < 2 the uniform models are all gap- 
less, and for the special values d = 2 cos they corre- 
spond to a sequence of conformally invariant Andrews- 
Baxter-Forrester (ABF) |1] models with central charges 
Cfe = 1 - 6/{k + l){k + 2) [9]. Physically, these ABF 
models can be thought of as describing chains of inter- 
acting non-Abelian particles described by su(2)fc Chern- 
Simons-Witten theory, believed to be relevant for certain 
quantum Hall states |10j . The known universal entangle- 
ment scaling of conformally invariant 1-fl dimensional 
systems [11] then implies that the entanglement entropy 
of a block of L contiguous sites, SJ^ , in the ground states 
of these models will scale logarithmically for L ^ 1 as 

When the Ji's are random, the Hamiltonians ([T|) can 
no longer be solved exactly. However, the RSRG ap- 
proach of Fisher [2] can be straightforwardly applied for 
all d > \/2 |5] with the result that the ground states 
all flow to the same inflnite randomness fixed point [T5] 
— one for which the bond strength distribution is the 
same as that of the fixed point of the random Heisenberg 
chain [2]. For this fixed point, Refael and Moore [13] 
have shown that if ul is the number of valence bonds 
leaving a given block of size L (see Fig. 1(c)), then, in 
the L > 1 Hmit, rT^ ~ ^^^^^ ~ ^ log^ L, where the 
overbar denotes a disorder average over random singlet 
states produced by the RSRG. This logarithmic scaling 
is a direct consequence of the inverse-square distribution 
of valence bond lengths characteristic of random singlet 
states [2]. Multiplying nZ by the entanglement per bond 
of log2 d then yields the RSRG result for the asymp- 
totic scaling of the entanglement entropy for the random 
ABF models, which is again logarithmic and has the form 
Sl^ ^rTZ log2 d ~ log2 L dlS]. 

In VBMC the ground states of the Hamiltonians ([ij are 
projected out by repeatedly applying ~H to a particular 



non-crossing valence-bond state \S). The result of this 
projection after n iterations is, 

i-Hr\s)= J2 J^,---A.K---K\s)- (2) 

The properties of the projection operators shown in Fig. 
1(b) imply that H^^ • --Il^JS) = Xi^^...^ija) where |a) is 
a non-crossing valence-bond state with the same norm 
as \S) and Xi^^...^i^ = d~™ where m is the number of 
times a projection operator acts on two sites which are 
not connected by a valence bond when projecting \S) 
onto |a). This projection thus leads to an expression for 
the ground state which becomes exact in the limit of 
large n (in our simulations we find it is sufficient to take 
n = 20N where TV is the number of sites) and has the 
form 

|V)-5]«.(«)|a), (3) 

a 

where w{a) — Ji^ - ■ ■ Ji^Xi^...i^. In VBMC the valence- 
bond states I a) contributing to |^) are sampled with 
probability w(a) by updating the sequence of projection 
operators (ii, • • • , i„) using the usual Metropolis method 

Given any observable O with expectation values 
0{a) — {a\0\a) / {a\a) in the non-crossing valence-bond 
states I a), VBMC can be used to compute the aver- 
age (O) = ^(q^)O(Q!)/ ^("^) fo'^ ^'^y state Ii/j) of 
the form ([3|, provided w{a) > 0. In what follows, an- 
gle brackets will always denote this average, though it 
should be noted that (O) will in general not be equal 
to (i/'|0|'(/')/(V'IV'); both because the valence-bond states 
are nonorthogonal and because the weight factors w(a) 
are amplitudes and not probabilities. One such averaged 
quantity is the valence-bond entanglement entropy, Sj^^ , 
which, for the uniform Heisenberg chain, is defined to be 
equal to the average number of valence bonds leav- 

ing a block of size L [HI [16] . To generalize to the 
ABF models with d = 2 cos it is natural to multiply 
('^l) by the asymptotic entanglement per bond of log2 d. 
For this choice, provided ^ 1, will be equal to 
S]^ for any single valence-bond state. We therefore take 
5r = (nL)log2 d. 

While is easy to compute numerically by VBMC, 
for a general superposition of valence-bond states it 
will not be equal to SJ^ . Nonetheless, VBMC simula- 
tions [ini dni dZ] of the uniform AFM Heisenberg chain 
with ~ 100 spins have shown numerically that 
grows logarithmically with L, in the same fashion as 
the von Neumann entanglement SJ^ . To characterize 
this log scaling it is convenient to introduce an effec- 
tive valence-bond central charge, c^^, defined so that 

1 ro VB 

Sl^ ~ ^ log2 L in the limit L > 1. 

In addition to showing log scaling of S'J^, previous 
VBMC simulations of the uniform AFM Heisenberg chain 
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FIG. 2: (Color online) Log-linear plots of valence-bond entan- 
glement as a function of block size L for uniform (upper 
panel) and random (lower panel) models with fc = 2, 3 and 
oo. For uniform models the solid lines correspond to the exact 
asymptotic scaling which follows from [18] . For random mod- 
els the solid lines show the asymptotic scaling predicted by 
the RSRG [13] . Results are for periodic chains with = 1024 
sites. For random chains the Ji's are taken to be uniformly 
distributed in the interval [0, 1] and results are self-averaged 
over all blocks for 50 disorder samples. 




FIG. 3: (Color online) d dependence of the valence-bond cen- 
tral charge, c^^, for uniform and random chains, along with 
the true central charges for the conformally invariant uniform 
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in units of logj d — the entangle- 



The solid red line is the exact result of Jacob 
[18| for {ul) (which equals c^^ 



I logj d when 



d = 2 cos j^) vs. d for uniform models, the black squares are 
the central charges of the ABF models for different values of 
k (the dashed line is a guide to the eye), the green line is the 
RSRG result c^^/log2d = ln2 for the random models, and 
the blue circles and triangles are our VBMC results for uni- 
form and random models, respectively. The inset shows the 
finite size extrapolation used to find c^^ for the uniform mod- 
els with fc — > oo and k — 2. Results were obtained by fitting 
S^-^ to ao + (c^-^/3) log2 X over the range N/16 < L < iV/4, 
where x — (N/tt) sm{Ln /N) is the conformal distance on the 
chain. The extrapolation clearly show the strong finite size 
effects for the k ^ oo case with d — 2. 



have given results consistent with c^^ being close to [16] , 
or even possibly equal to [15 , 1, the value of the true 
central charge for the uniform d — 2 model. However, 
Jacobsen and Saleur [T5] were able to determine the exact 
asymptotic scaling of (n^) analytically for all d <2, and 
while their results confirmed the log scaling of S^^ for 
L 3> 1, they also showed that for d = 2 the coefficient 
for this scaling does not correspond to c^^ being equal 
to one, but rather c^^ = 121n2/7r2 ~ .843. 

Figure llJa) shows our VBMC results for Sj^ for 
k = 2,3 and oo (corresponding to d = \/2, <P and 2, 
respectively, where (f> is the golden mean) for periodic 
systems with N — 1024 sites. The solid lines show the 
exact asymptotic scaling found by Jacobsen and Saleur 
[18j which clearly agree with our numerical results for 
L ^ 1. Note that for the case fc cx) it is necessary 
to consider fairly large values of L before entering the 
scaling regime, whereas for k = 2 and 3 the scaling be- 
gins at relatively small L. This fact may account for the 
initial numerical difficulty in determining c^^ for d = 2 
using small systems (see, however, |17j). Presumably, 
the reason finite size effects become more pronounced as 
d approaches 2 is because this is a critical value (for d > 2 
the uniform models acquire a gap [7] ) . 



For random Heisenberg chains was first computed 
numerically by Alet et al. [TS]. Following the same pro- 
cedure as these authors, we compute by determining 
{ul) for particular realizations of disorder and then dis- 
order averaging. For the random ABF models we again 
multiply by the entanglement per bond, and thus take 
= {til) log2 d (here again the overbar denotes dis- 
order average) . Figure |2jb) shows log plots of our results 
for SY^ for random chains, again for A: = 2, 3 and oo and 
N — 1024. The solid lines show the scaling predictions of 
Refael and Moore [T3] for based on the RSRG which 
clearly agree with our numerical results. As pointed out 
by Alet et al. [15], the fact that SY^ and Sl^ show the 
same scaling for L » 1 is to be expected if, as predicted 
by the RSRG, on long length scales the ground states of 
the random models are dominated by a single valence- 
bond state. 

Figure [3] shows our VBMC results for c^^ for both 
uniform and random models for different values of d. For 
the uniform models Fig. [3] also shows the exact values of 
c^^ which follow from the analytic results of Jacobsen 
and Saleur 18J as well as the true central charges Ck of 
the ABF models with d = 2 cos [TH]. For the random 
models the d dependence of c^^ is seen to be entirely due 
to the entanglement per bond, reflecting the fact that the 
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FIG. 4: (Color online) Log-linear plots of bond fluctuations as 
a function of block size L for both uniform and random chains 
with — > oo (upper panel) and fc = 2 (lower panel). For 
uniform chains these fluctuations grow logarithmically with 
L in agreement with the analytic results of 18 (solid lines), 
indicating that bonds are strongly fluctuating on all length 
scales. For random chains the fluctuations saturate, signaling 
the formation of a random singlet phase in which the bonds 
have locked into a particular random singlet conflguration on 
long length scales. The parameter u (deflned in the text) 
is a measure of the disorder strength and the results clearly 
show that the saturation length scale grows with decreasing 
u. Results are for periodic chains with N — 1024 sites and 
for random models are self-averaged over all blocks for 100 
disorder samples. 



valence bond length distribution, which determines the 
coefficient in front of the log scaling of (n^), and which 
depends on d for the uniform case, becomes independent 
of d when disorder is included. 

The expectation that for random chains the L 3> 1 
scaling of should be the same as that of S}^ is 
based on the assumption that the valence bonds lock into 
a particular random singlet state on long length scales. 
This assumption is in turn based on the RSRG approach 
which, although it can be shown to capture the long dis- 
tance properties of the fixed point exactly [2 , is still an 
approximate method. It is clearly desirable to have a di- 
rect numerical demonstration that the valence bonds are 
indeed locking into a particular random singlet state on 
long length scales. 

To provide such a demonstration we calculate the fluc- 
tuations in n^. To be precise, we first compute the 
quantity (n|) — (n^)^ for a particular block of size L 
and a particular realization of disorder, and then per- 



form a disorder average. The quantity we compute is 
thus (t£ = (ri|) — {ul)'^. For this choice of averaging 
aj^ has the property that, in an idealized random singlet 
phase for which the ground state is a single non-crossing 
valence-bond state, cr£ would vanish, even though the 
number of bonds leaving a given block would be different 
for different realizations of disorder. 

For the uniform models with d < 2 Jacobsen and 
Saleur [18] have also determined the asymptotic scaling 
of cr£ (in this case there is, of course, no disorder aver- 
age). They find that, like (ul), a\ scales logarithmically 
with L, and they were again able to obtain exact analytic 
results for the coefficients of the log. This log scaling 
characterizes the expected strong fiuctuations of valence 
bonds on all length scales in the uniform models. 

For random models there will be strong valence bond 
fiuctuations on length scales up to a certain length I, be- 
yond which the bonds should lock into a random singlet 
phase. In this case one expects a\ to not differ much 
from its value for the uniform case when L <^ I, but for 
L I the fluctuations should saturate. This saturation 
is due to the fact that, once the block size L becomes 
much larger than the fluctuation length I, the fluctua- 
tions which occur outside of a distance I from the two 
boundaries of the block will not change the number of 
bonds leaving the block, and hence will not contribute to 

Figure [4] shows log plots of our results for cr^ for the 
case k oo (corresponding to the Heisenberg chain, with 
d = 2) and fc = 2 (corresponding to the Q — 2 Potts 
model, i.e. the critical transverse field Ising model, with 
d = y/2). For the uniform models our results confirm the 
log scaling and exact coefficients found by Jacobsen and 
Saleur [TB]. For the random models the Ji's are taken 
to be uniformly distributed in the interval [1 — w, 1 -f u], 
where m is a measure of disorder strength, and we ob- 
serve that (t£ saturates on a length scale which grows as 
u decreases. We believe the observation of this satura- 
tion, which indicates a finite fiuctuation length scale I, 
together with the log scaling of (ni), which indicates a 
power-law distribution of valence bond lengths, provides 
a direct numerical proof of random singlet phase forma- 
tion in these models. 
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